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I. INTRODUCTION 

Is there any new physics (NP) beyond the Standard 
Model (SM)? Probably yes, if high-energy physics does 
not want itself to be lost in a mire of unanswered ques- 
tions. Supersymmetry (SUSY) is one of the most widely 
discussed options of NP, in both its R-parity conserv- 
ing (RPC) and R-parity violating (RPV) incarnations. 
However, SUSY introduces a plethora of new particles, 
and even in its most constrained version, a few more arbi- 
trary input parameters over and above to that of the SM. 
Thus, it has become imperative to constrain the SUSY 
parameter space as far as possible from existing data. 

There are, of course, direct bounds on the sparticle 
masses from collider data. These bounds are typically 
weaker for RPV SUSY than for the RPC version, since 
in the former case the final-state signal is radically dif- 
ferent and more possible channels are open. Apart from 
them, there are bounds on the parameter space from low- 
energy processes, whose amplitudes can be affected by 
intermediate sparticle states. For RPC SUSY, the new 
amplitudes must be at least at one-loop level, so that they 
can compete with the SM amplitude and show up only 
if the SM amplitude is also suppressed. The new am- 
plitudes may appear at tree-level if R-parity is violated, 
and hence there is a greater chance of constraining RPV 
SUSY models from low-energy data. There is another 
good reason to focus upon the RPV version: most of the 
precision low-energy processes from which one can obtain 
bounds are flavor-changing neutral current (FCNC) type. 
Now, there are some elegant mechanisms of suppress- 
ing FCNC effects in RPC SUSY, but absolutely none 
(except putting some flavor-changing Yukawa-type cou- 
plings equal to zero, or vanishingly small, by hand) if we 
have RPV. Thus, the bounds on the parameter space of 
RPV SUSY coming from FCNC processes are in a sense 
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more robust. 

In this paper we use the data from — K'^ and 
— B'^ mixing to constrain the relevant couplings for 
RPV SUSY. For the latter we use both Ams and sin(2/3) 
constraints, while for the former we use the results on 
ArriK and Sk- We do not discuss other CP violating pa- 
rameters like e' /e, since that has large theoretical uncer- 
tainty. We have also discussed some correlated channels 
leading to leptonic and semileptonic decays. We do not 
consider the Bg system since there is only a lower bound 
on AMb, ■ This, in turn, means that there is no such 
upper bound on the relevant NP couplings; we can only 
have a lower bound, which is consistent with zero. The 
situation should change dramatically once the hadronic 
B machines, producing copious Bg mesons, come on line. 

Do we have any motivation to invoke NP for the K 
and the B systems? In other words, is there any incon- 
sistency of the experimental data with the SM predic- 
tions? The answer is yes for the B system, though the 
error bars are still large to draw any definite conclusion 
(but we have reasons to be hopeful). The sore thumbs 
are (i) the abnormally high branching ratios (BR) for 
the generic channels B rj'K,r}K* [1] [32], (ii) the di- 
rect CP-asymmetry in the channel Bd tt'^tt^ as found 
by Belle [2], (iii) the discrepancy in the extracted value 
of sin(2/3) from Bd — > J/4'Ks and Bd — ^ (f'Ks [3], and 
(iv) the so-called "ttK" puzzle: the abnormal enhance- 
ment of electroweak penguins in B ttK decays [1] 
[33]. However, one must not be over-enthusiastic since 
these channels are nonleptonic and QCD uncertainties 
are yet to be fully understood. But one may hope more 
such anomalies from leptonic and hadronic B-factories. 
For the K system, the nonleptonic channels are notori- 
ously difficult for any systematic analysis of NP effects 
[4], but for the first time we are having precise data (or 
bound) on leptonic and semileptonic K decay channels 
from Brookhaven and Da$ne. It is always better to be 
ready for any unexpected result. 

Another motivation for supersymmetry comes from the 
neutrino mass. It has been shown that the existence of 
nonzero neutrino mass may have observable signatures in 
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B-physics [5, 6] like the enhancement of the b sss am- 
phtude and also AMg, the mass difference between two 
Bs mass eigenstates, if one assumes some unified RPC 
SUSY theory. In RPV SUSY, one may have a definite 
texture at the GUT scale (only a few RPV couplings 
are nonzero) and can generate the whole lot of nonzero 
RPV couplings at the weak scale through renormalization 
group evolution and CKM-type mixing [7]. One expects 
a signature of those RPV couplings relevant for neutrino 
mass generation in the B-factory data. 

Must the RPV couplings be complex? The answer is 
that there is no a priori reason why they should be all 
real. Even if there is some real GUT texture, the weak 
scale couplings may turn out to be complex. The phase 
of one single coupling can be absorbed in the sfermion 
field, but a nontrivial phase should be there in a product 
of two such couplings. 

Effects of RPV SUSY on K and B physics have been 
discussed extensively in the literature [8-12]. Constraints 
coming from K° — and B*^ — mixing have been 
discussed in [13-15]. This paper is a culmination of the 
effort in the following sense: (i) This is the first time 
that both SM and all RPV contributions (including their 
phases) have been taken into account for the systems 
under consideration; (ii) the short-distance QCD correc- 
tions have been implemented upto next-to-leading order 
(NLO) accuracy (there is scope for questioning this pro- 
cedure for the K system; this will be discussed later); (iii) 
Both CP-conserving and CP-violating constraints have 
been discussed for the B'^ — B° and — systems; 
(iv) Some discrepancies in the relevant formulae used in 
some of the earlier works have been corrected; and (v) 
Correlated signals for leptonic and semileptonic decays 
have been discussed. 

The paper is arranged as follows. In Section 2 we out- 
line the relevant formulae necessary for the analysis. Sec- 
tion 3 deals with the numerical inputs and Section 4 with 
the analysis of the results. Decay channels mediated by 
the same couplings are touched upon in Section 5, while 
we conclude and summarize in Section 6. Some calcula- 
tional details have been relegated to the two appendices. 



II. BASIC INPUTS 
A. Neutral meson mixing 

L et th e neutral mesons be generically denoted by M° 
and M^, with the valence quark content qd and qd re- 
spectively. For the cases under study, q can be either d 
or s. For Bg system, replace qhy b and dhy s. 

The off-diagonal element in the 2x2 effective Hamil- 
tonian causes the Af° — mixing. The mass difference 
between the two mass eigenstates AM is given by (fol- 
lowing the convention of [16]) 



with the approximation |Mi2| ^ |ri2|. This, however, is 
true for the B system only. Let the SM amplitude be 



Mf2^|exp(-2ie5M) 



(2) 



where Osm = f3{4>i) for the — system and approx- 
imately zero for the — (and also for Bg — Bg) sys- 
tem. We follow the (a, (3, 7) convention for the unitarity 
triangle [16]. 

If we have n number of NP amplitudes with weak 
phases 6n, one can write 

n 

M12 = \Mff\ exp{-2iesM) + J2 IMlsl exp(-2ieO- (3) 

4 = 1 

This immediately gives the effective mixing phase 9eff 
as 

_1 [Aff^^l sm{20sM) + Ei \Ml2\ Mm 

2 ^'^'^^ [Mf^^l cos(2^sm) + \Mh\ cos(20,) ' 

(4) 

and the mass difference between mass eigenstates as 
AM = 2[|Mf2^|2-j-^|M^2|2 

i 

+ 2|Mf2^ I 1 cos 2{esM - ei) 

i 

+ 2^^|Mf2||M^2|cos2(0,-e,)]'/'- (5) 

i j>i 

These are going to be our basic formulae. The only task 
is to find MI2 and Oi. 

For the - K° system [17], |ri2| is non-negligible, 
and we can write 



AM 
AT 



2 Re 
-4 Im 



(M12 - |ri2)(Mi*2 - |r^2) 



1/2 



(M12 - ^ri2)(Mr2 - ^ri2) 



1/2 



,(6) 



so that AM = — (l/2)Ar. Since the dominant decay is 
to the / = final state, one can neglect /mri2 and write 



AM = 2 Re M12, Ar = 2 i?e ri2. 
The CP-violating parameter ek is given by 

1 Im M12 



which can be written as 



2^/2 Re M12 ' 



1 Im M12 
71 AM • 



(7) 



(8) 



(9) 



AM = 2|Mi 



(1) 



Note that Re M12 has both short-distance (SD) and long- 
distance (LD) contributions. The LD contribution is not 
calculable; what one calculates from the box amplitude 
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is the SD part. That is why one generally uses the ex- 
perimental value of AM/f in the denominator of eq. (9). 
For — system, the short-distance SM amplitude 



IS 



(i^O|i/e//|i^°) 



2m 



K 



^{VcciV:^Ym'mKfKBKmi^So{x,),{lQ) 

where generically Xj = m^^ /ra^^, fx is the K meson decay 
constant, and rjx (also called r]cc in the literature) and 
Bk parametrize the short- and the long-distance QCD 
corrections, respectively {B^ is a phenomcnological pa- 
rameter, thrown in to parametrize the LD contribution). 
The tiny top-quark loop dependent part responsible for 
CP violation has been neglected. The function is given 
by 



^o(^) = 



4x - llx'^ + x^ 3x^ In X 



4(1 -xf 



2{l-xf 



(11) 



For the S° — 5" system, we have an analogous equation, 
dominated by the top quark loop: 

SM _ mHeff\B<^) 



2mi 



^{VtdVtlYriBmBfBBBmi^Soixt). (12) 

There is not enough motivation to consider the Bs — Bg 
system right now, since there exists only a lower bound on 
AMs(> 14.4ps~^) [1]. This can accomodate arbitrarily 
large NP couplings, starting from zero. However, if the 
NP amplitude has a nonzero phase, then there will be an 
effective phase in Eg — Bs mixing amplitude, whose pres- 
ence may be tested in the hadronic B factories. Whether 
there is any detectable new physics in Bs — Bs mixing 
(this is particularly relevant since the b s penguin 
transition, at least for nonleptonic decays, shows hint 
of an anomalous behaviour) can be effectively tested in 
hadronic B machines, which will measure AMs as well 
as the CP asymmetries in Bs J/tpij^ and hopefully 
Bg — > (t){r]')(f){r]'). Detection of NP signals in the double 
vector meson modes rcqiiires angular analysis of the de- 
cay products but hopefully can be done in future colliders 
[18]. 

In the presence of NP, the general AF = 2 effective 
Hamiltonian can be written as 



-i/AF=2 
^eff 



where /j. is the regularization scale, and 

02 = {qPRd)i{qPRd)i, 

03 = iqPRd)siqPRd)s, 
Oi = {qPLd)i{qPRd)u 
O5 = {qPLd)s{qPRd)s, 



where q is either b or ,s, and Pr(l) = (1 + (^)75)/2- The 
subscripts 1 and 8 indicate whether the currents are in 
color-singlet or in color-octet combination. The ds are 
obtained from corresponding OjS by replacing L ^ R. 
The Wilson coefficients Cj at q"^ = include NP ef- 
fects, coming from couplings and internal propagators. 
However, for most of the NP models, and certainly for 
the case we are discussing here, all NP particles are heav- 
ier than mw and hence the running of the coefficients 
between mw and /U = 0{ms or m^) are controlled by 
the SM Hamiltonian alone. In other words, NP deter- 
mines only the boundary conditions of the rcnormaliza- 
tion group (RG) equations. For the evolution of these 
coefficients down to the low-energy scale, we follow Ref. 
[19], which uses, for B° — B^ mixing, n = mb = 4.6 GeV. 
The low-scale Wilson coefficients, using the NLO-QCD 
corrections, are 

= E E i^r' + ^^^'1 V'^'CsiMs) (15) 



where r] = as{Ms)/as{mt), Ms being the scale of NP, 
which, for SUSY, may be taken to be the average of the 
squark and the gluino masses. We will use Ms = 500 
GeV throughout the paper. For the numerical values 
of a, 6 and d matrices we refer the reader to eq. (10) 
of Ref. [19] (to avoid confusion with the Wilson coeffi- 
cients, we use the symbol d, for the matrix denoted by 
c in [19]). All the numbers are not relevant for our dis- 
cussion; the modification of the NP operators due to the 
short-distance QCD corrections is discussed in Section 
4. It will be shown that only operators Oi and O4 are 
relevant at the scale q"^ = m^, while there is a slight 
admixture of O5 at the low-energy scale. 

The operators Oi are also to be renormalized at the 
scale fi. The expectation values of these operators be- 
tween B° and B^ at the scale /j. are given by 



^Ci{iJ,)Oi{iJ,) + ^Ci{ii)di{iJ,) + H.c. (13) where 



(Oi(m)) 

(02(/i)) 
(03(m)) 

(05(m)> 

Sb = 



-mBfBBiin), 

-^SBmlfMli), 

^SBTnlflBiiii), 
7:SBm\f%B^{lj), 



n^B 

mb{mb) + mdimb) 



(16) 



(17) 



(14) 



The B-parameters, whose numerical values arc given in 
Section 3, have been taken from [20]. Note that the ex- 
pectation values are scaled by factor of 2mB over those 
given in some literature due to our different normaliza- 
tion of the meson wavefunctions. It is trivial to check that 
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both conventions yield the same values for physical ob- 
servables. Analogous expressions follow for the K° — K° 
system, with mB^fs triK^fK, and 

S^^-( y, (18) 

with n = 2 GeV. The running of the Wilson coefficients 
is given by an equation which is exactly analogous to eq. 
(15) [21] but with different a, b and d matrices, which 
are evaluated at ;U = 2 GeV. The -B-parameters, which 
include all nonperturbative effects below the scale /U, are 
different too. 

It is noteworthy that the NLO corrections as discussed 
above cannot be applied when there are light quarks like 
u flowing in the loop, since the operator product expan- 
sion takes a different form [21]. This will be the case 
when we switch on the RPV interaction. There should 
be an enhancement coming from large logarithms; how- 
ever, this puts on a tighter constraint on the RPV cou- 
plings. Since these QCD corrections are not known, we 
use the same procedure as adopted for heavy particles, 
and probably deduce somewhat lenient constraints on the 
parameter space for those particular couplings. 

It is clear from equations (10) and (12) that at least in 
the SM, one can either evaluate the matrix element of Oi 
at niw and compute the running to the low energy scale, 
or alternatively can use the parameters Bk,b and r]K.B, 
to be determined from the lattice studies. We adopt the 
latter method for calculating the SM amplitudes, scan- 
ning over the whole range of those parameters, so that 
even the nonperturbative effects can be successfully taken 
into account, including all the uncertainties. This gives 
us the most conservative bounds on the NP parameters. 
The NP amplitude is calculated using the former method. 

B. R-Parity violating SUSY 

R-parity is a global quantum number, defined as 
(_]^-)3B+i+2S^ which is -1-1 for all particles and -1 for 
all superparticlcs. In the minimal version of supcrsym- 
metry and some of its variants, R-parity is assumed to 
be conserved ad hoc, which prevents single creation or 



With such a term, one can have two differen t kin d of 
boxes, shown in Fig. 1 , that contribute to M° — mix- 
ing: first, the one where one has two sfermions flowing 
inside the loop, alongwith two SM fermions [14], and scc;- 
ondly, the one where one slepton, one W (or charged 
Higgs or Goldstone) and two up-type quarks complete 



annihilation of superparticlcs. Such models have a nat- 
ural candidate for dark matter, and have telltale signa- 
tures in collider searches as large missing energy. How- 
ever, models with broken R-parity can be constructed 
naturally, and such models have a number of interesting 
phenomenological consequences [22, 23]. Some of these 
R-parity violating models can be motivated from an un- 
derlying GUT framework [24]. 

It is well known that in order to avoid rapid proton 
decay one cannot have both lepton number and baryon 
number violating RPV couplings, and we shall work with 
a lepton number violating model. This leads to both slep- 
ton (charged and neutral) and squark mediated decays, 
and new amplitudes for — mixing, where parti- 
cles flowing inside the box can be (i) charged slcptons and 
up-type quarks, (ii) sneutrino and down type quarks, (iii) 
squarks and leptons. One or both of the scalar particles 
inside the box can be replaced by W bosons, charged 
Higgs bosons and Goldstone bosons (in a non-unitary 
gauge) (see Fig. 1). We follow the usual practice of avoid- 
ing the so-called "pTirc SUSY" contributions to the box 
amplitudes, i.e., those coming from charginos, ncutrali- 
nos or gluinos inside the loop. Not only the strongly in- 
teracting superparticlcs are expected to be heavier than 
the electroweak ones (and hence the contribution be- 
ing suppressed), but also one can choose SUSY models 
where these contributions become negligible {e.g., align- 
ment in the squark sector, or Higgsino-dominated lighter 
chargino, to kill off the respective boxes.) Since the cur- 
rent lower bound on the slepton mass is generally weaker 
than that on squark mass by a factor 2-3, the slepton 
mediated boxes have greater chance to be numerically 
significant. 

We start with the superpotential 

Wv = KjkLiQjDl (19) 

where i,j,k = 1,2,3 are quark and lepton generation 
indices; L and Q are the S'f/(2)-doublet lepton and quark 
superfields and is the S'L''(2)-singlet down-type quark 
suporfiold respectively. Written in terms of component 
fields, this superpotential generates six terms, plus their 
hermitian conjugates: 



+ H.c. (20) 
I 

the loop [13]. It is obvious that the first amplitude is 
proportional to the product of four A' type couplings, 
and the second to the product of two A' type couplings 
times Gf- We call them L4 and L2 boxes, respectively, 
for brevity, where L is a shorthand for A'. 

We will constrain only products of two A'-type cou- 
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plings at a time, and assume a hierarchical structure, 
i.e., only one product is, for all practical purpose, si- 
multaneously nonzero (but can have a nontrivial phase) . 
This may not be physically the most appealing scenario 
but keeps the discussion free from unnecessary compli- 
cations. For any product, there are individual bounds 
on each of them, the product of which we call the di- 
rect product bound (DPB). Interesting bounds are those 
which arc numerically smaller, and hence stronger, than 
the corresponding DPBs, the more the better. This turns 
out to be the case for almost all the products. A list of 
the individual bounds at the weak scale can be found in 
[7]. 

At this point, let us clarify how the bounds are ob- 
tained. There are, generically, three boxes: SM, L2 and 
L4. For L2 boxes, one can have same or different up- type 
quarks flowing inside the loop. For the former case, there 
is also an L4 counterpart, i.e., the coupling that gives rise 
to the L2 box can also generate an L4 box. The operator 
for the L2 box is O4, while that for the L4 box is Oi. 
The RG evolution generates a small O5 admixture to O4 
at the low-energy scale, while O1 is just multiplicatively 
renormalizcd. 

The effective Hamiltonian is obtained from a coherent 
sum of these three amplitudes, and there can be intri- 
cate interference patterns depending on the magnitudes 
and the phases of them. A general trend, however, is 
easy to follow. The pure SM part is proportional to G|, 
times some CKM factor. The L2 box is proportional to 
GpX^ jvr?, and the L4 box goes as /m^. For small 
values of the A' A' coupling, L2 contribution dominates 
over L4, and the "bound" (as obtained by [13]) is con- 
trolled by L2. We will see that if the phases are included, 
this procedure docs not give the bounds, so the munbcrs 
obtained in this way are not listed in tables 1 and 2. For 
larger values of the coupling, L4 becomes dominant, and 
if there happens to be a cancellation between SM and 
L2 amplitudes, L4 controls the show. That is, the solu- 
tion describes a contour in the Re{X'X') — Im{X'X') plane. 
There are certain points, far from the origin, where there 
happens to be a complete destructive interference be- 
tween SM and L2 boxes. Thus, if the product has an 
arbitray phase, it is always the L4 term that gives the 
bound when the quarks inside are the same. We will see 
explicit examples of this later. 

If the up- type quarks in the box are different, there is 
no corresponding L4 contribution (actually, there is, but 
that comes from four different A' type couplings which 
we neglect systematically), and the treatment becomes 
simpler. 

A note of caution here. In the literature, most of the 
bounds on the RPV couplings are not absolute, but scale 
with the corresponding sfermion mass that comes in the 
propagator. For a product coupling, the dependence is 
typically A'A'/wj. The results are usually quoted for 
= 100 GeV, and bounds for higher are straight- 
forward to obtain. However, for the L4 boxes, there can 
be two distinct type of amplitudes, one through a quark- 



slcpton box, and the other through a squark-lepton box. 
For our analysis, we take all sleptons to be degenerate at 
100 GeV, and all squarks to be degenerate at 300 GeV 
(100 GeV squarks are already ruled out). Thus, it is 
not easy to scale the bounds quickly; a good way is to 
remember (and scale accordingly) that at the amplitude 
level, the squark contribution is about 10% of the slepton 
contribution, and the second way, to get a more rough 
scaling, is to neglect the squark contribution altogether. 

The effective Hamiltonian for the L4 box has been com- 
puted in [14], with degenerate sfermions. For a nonde- 
generate case, the expression reads {q = 2 for K° — 
and g = 3 for 5° - 5° boxes) 



Wi4 = 



iKklKkqf 

1287r2 




where 



I{x) = 



1 — + 2x log x 
{l-xf 



Oi 



(21) 



(22) 



with /(O) = 1, which justifies the omission of this func- 
tion in the last term, where we have explicitly assumed all 
leptons to be massless. In fact, this is also a reasonable 
assumption for all quarks except the top quark. 

There may be another L4-type amplitude, with left- 
handed quarks as external legs and right-handed down- 
type squarks and neutrinos (and quarks and sneutri- 
nos too) flowing inside the box. Neglecting the fermion 
masses, the effective Hamiltonian reads 



1287r2 



1 



1 

ml 



Oi 



(23) 



Note that this bounds a different combination. As we 
will see soon, there is a corresponding L2 box with the 
same coupling as in Hl4, but none for ^^4. (The reason 
is simple: VK-bosons couple only to Icft-chiral quarks, 
and scalar couplings are negligible.) Thus, all A^^^A^^j. 
combinations are expected to yield the same bound. 

The case for the L2 box is more complicated. A sim- 
plified treatment, keeping all bosons (sleptons, Ws, and 
Higgs bosons) to be degenerate at 100 GeV, neglecting 
the QCD corrections, and saturating the experimental 
number with RPV amplitude alone, was given in [13]. 
Only real RPV couplings were considered since there 
was no data for CP violation in the B system at that 
time. We improve the calculation by incorporating all 
the above factors, but note that the bounds quoted in [13] 
for tan (3 = 1 are still surprisingly close to the mark. In 
fact, with the introduction of RPV phases, the possibility 
of destructive interference between the amplitudes opens 
up, and one would have expected a weaker bound than in 
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FIG. 1: R-parity violating contributions to B° — B° mixing. Figure (a) corresponds to L4, while figure (b) to L2 amplitudes 
(see text for their meanings). For L4, there are similar diagrams with squarks and leptons (both charged and neutral), as well 
as diagrams with left-chiral quarks as external legs and quarks and sneutrinos flowing in the box. For L2, there are diagrams 
where the W is re-placed by the charged Higgs or the charged Goldstone. The internal slepton can be of any generation, and so 
can be the internal charge +2/3 quarks, generically depicted as u. 



[13]. Consideration of the CP violating parameters sin 2(3 
and £k helps us to put a tighter set of constraints. The 
effective Hamiltonian, in this approximation, is given by 

_ Gp^ iki^' ikq ^r* \r s, 

n,. - — ^^f^^k^VM X 

[(1 + cot^ l3)xlJ{xk) + I{xk)\ O4, (24) 
where Xk = rn^^/rn^, and 

-2{x -l) + {x + l)\ogx 



J{x) = 



{x-lf 



(25) 



and cot/3 = Vd/vu, the ratio of the vacuum expecta- 
tion values of the two Higgs bosons that give mass to 
the down- and the up-type quarks respectively (not to 
be confused with the phase of Vtd)- For nondegenerate 
masses the expression is expectedly more complicated: 



Gp^'ikl^'ikg 



(26) 



where ^^^'*^s are factors that come out of the box di- 
agram integration (the superscripts indicate the SM bo- 
son in the loop). Note that as expected, the Higgs- and 
Goldstone-box contributions are negligible for k = 1,2. 
The expressions for these factors in 't Hooft-Feynman 
gauge are given in Appendix A. 

The Hamiltonian is slightly modified if we consider two 
different up- type quarks k and p inside the loop: 



'Hl2 



Gp^'iki^'i 



4\/27r2 



lElV* V J •>< 



[m2j,Sf2 + cot^ PmlmlB^^ + mlmlB%] O4, 



(27) 



and again the integration factors 3^2^'^ are shown ex- 
plicitly in Appendix A. Note that the Hamiltonian is 
evaluated at q"^ = and the matrix elements are to 
be evaluated at the proper low-energy scale. 

Apart from these new amplitudes, there; are a num- 
ber of other boxes, a detailed list of which is given in 
Appendix B , which arc proportional to four A' type cou- 
plings. They may be important in a scheme with a defi- 
nite texture for such couplings, but for our present study, 
we do not consider them any further. 



III. NUMERICAL INPUTS 

The major sources of the numerical inputs are: (i) the 
Heavy Flavor Averaging Group (HFAG) website [1] for 
the latest (summer 2003) updates on B physics; (ii) Par- 
ticle Data Group 2002 edition [25] and update for 2004 
available on the web [26] ; and (iii) the inputs used in the 
CKMfittcr package [27]. The quark masses and Wilson 
coefficients have been taken from [19, 21]. We use the 
following numbers. 

Masses (all in GeV) [19, 21, 27]: 

ms = 5.2794; mx = 0.494; rUb = 4.23; 

mt = 167 ±5; rUc = 1.3; ms(2 GeV) = 0.125; 

rud = 0.007. (28) 

The quark masses have been evaluated in the MS 
scheme. The pole mass for the top quark is about 5 GeV 
higher and the mass for the bottom quark is 4.6 GeV. 
The strange quark mass, however, is a major source of 
error for the chiral ly e nhanced operators O4 and O5 eval- 
uated for the K^ — K'^ system. The function Sk increases 
by a factor of 2 if we take nis = 95 MeV, which in turn 
tightens the bound on the A' A' couplings, nig = 125 MeV 
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is almost at the uppermost limit of the allowed range and 
hence generates the most conservative bound. 
The mass differences, in ps~^, are given by [1, 25] 

AruB = 0.502 ± 0.006; Atuk = (5.31 ± 0.01) x 10"^. 

(29) 

The CP-violating parameter sin(2/3) is taken as the 
average over all charmonium modes [1], since we have 
reasons to suspect that there may be hints of new physics 
in 6 ^ s transition. On the other hand, the theoretical 
prediction of sin(2/3), obtained from a fit excluding the 
direct experimental results, is taken from [28]: 

sin(2/3)exp = 0.736 ±0.049; 
sin(2/3)tfc = 0.685 ±0.052. (30) 

We also use [27] 

IskI = (2.282 ± 0.017) x 10"^. (31) 

This brings us to the only possible caveat of our anal- 
ysis. If we assume the existence of NP in the — 
mixing amplitude, one should not use Atob data for ex- 
tracting Vtd- In fact, what is needed is a rcevaluation 
of the bounds excluding all inputs that may be affected 
from NP {e.g., Athb, ^k, etc.). We do not venture into 
that here. Rather, we take the 90% confidence limit (CL) 
bound on Vtd assuming unitarity of the CKM matrix and 
considering tree-level processes only [26] : 

|14^| = 0.0094 ±0.0046. (32) 

Rest of the CKM elements are taken from [26], most of 
them at their central value: 

\Vud\ = 0.9745, \Vus\ = 0.224, IK^I = 0.0037(8), 

\Vcd\ = 0.224, \Vcs\ = 0.9737, \Vcb\ = 0.0415, 

\Vts\ = 0.040(3), \Vtb\ = 0.99913. (33) 

The angle 7 is taken to lie between 50° and 72° [27]. The 
analysis is fairly insensitive to its precise value. 

The long- and short-distance QCD corrections to the 
box amplitudes are mostly determined from lattice stud- 
ies. We use [27] 

Bk = 0.86 ±0.14 ±0.06, 
s/b^Jb = (0.228 ± 0.033) GeV, 
fK = 159.8 MeV, 

rtK = 1.38 ± 0.53; 773 = 0.55. (34) 

The leptonic and semileptonic BRs for the K meson, 
which are of interest to us, are as follows [26]: 

Br(if+ ^ n+vV) = 1.6j^-^ x 10"^°; 

Bt{K+ n+ ^l+ = (8.1 ± 1.4) X 10"*; 

Br(X+ ^ TT+e+e") = (2.88 ± 0.13) x 10"^; 

Br(ii:L ^ /X+/X-) = (7.24 ± 0.14) X 10'^ (35) 



There are stringent bounds on e~^e~ too but that 

is chirally suppressed compared to the fJ-^ IJ-^ mode. For 
the B mesons, the relevant numbers may be found in 
[1]. We will not analyze them here (see [29] for such 
an analysis) but the bounds on RPV couplings obtained 
from B decays are not compatible with those obtained 
here, particularly for a 300 GeV squark. 

To evaluate the QCD corrections, we take as{m%) = 
0.1172 ± 0.0020 [25], and take the SUSY scale Ms = 500 
GeV. The precise value of this scale is not important, 
however, and we can take it to be at the squark mass 
scale (300 GcV) without affecting the final results. The 
exact evolution matrix can be found in [21] and [19]; for 
our purpose, it is sufficient to note that for the K sys- 
tem, the operator Oi is multiplicatively renormalized by 
a factor 0.794 at the scale ji = 2 GeV, and the opera- 
tor O4 at mw changes to (3.96504 + O.I49O5). For the 
B system, the respective numbers are 0\ 0.820Oi, 
O4 (2.83O4 ± O.O77O5). We again stress that the- 
oretically the procedure is questionable for boxes with 
light quarks flowing in the loop. However, the numbers 
that we obtain are fairly robust and one can very well 
drop the NLO corrections altogether, if necessary, with- 
out compromising the results. 

The relevant B-parameters (eq. (16)) are [20, 21] 

For B : 

Bi{mb) = O.S7{4)tl, 
Bi{mb) = 1.16(3)tf, 
Br,{mb) = 1.91(4)tf , 

For K : 

Bi(Ai) =0.60(6), 
Bi{ij) = 1.03(6), 

B5(/i) = 0.73(10). (36) 

Since the O5 admixture is small, one can take the cen- 
tral values for these parameters without introducing too 
much error. 

In the supcrsymmetry sector, wc assume only one RPV 
product coupling to be nonzero at a time. We take all 
sleptons to be degenerate at 100 GeV, and all squarks 
at 300 GeV. We also take tan/3(= U2/W1) = 5 (very low 
values are excluded by LEP, and the numbers are not 
sensitive to the precise choice of tan /3) , and the charged 
Higgs boson mass as 200 GeV (lower values are disfavored 
from h S7). 



IV. ANALYSIS 

Our bounds are summarized in tables 1 and 2, for the 

K and the B systems respectively. The procedure has 
been outlined in Section 2, and is identical for both sys- 
tems. However, for the — system, the SD contri- 
bution may not be the full story. If the SD contribution 
is smaller, the contribution to Auik is also smaller, and 
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TABLE I: Bounds on A'A' combinations from K°-KO mrxme 
and ek, shown for cases when the product is strictly real and 
when the product may be complex. The table displays the 
magnitudes only, and not the signs. Weakest DPBs, displayed 

in the last column, occur for i = 3. They arc from [7] except 
that we scale them for 300 GeV squarks wherever necessary, 
unless they hit the perturbative bound. The entries marked 
with a dagger are bounded from ArriK, see text. The one 
mar ked with a double dagger is bounded from tree-level K" — 
K" mixing, and those marked with a 4|k have bounds of the 
order of 10"^ from squark-mediated —> -k^uV (see text). 



hence the analysis runs for a smaller Amx and gives 
tighter bounds on the relevant RPV couplings. The 
bounds we obtain are, therefore, the most lenient ones. 

It appears that the constraints on the — system 
arc essentially controlled by ek-, particularly when we 
allow for complex RPV couplings. In the SM, the CP vi- 
olating part is suppressed by the CKM factors; with RPV 
there is no such suppression, and thus one gets tight con- 
straints. They are shown in the third and fourth columns 
of table 1. The second column shows the constraints for 
real RPV couplings. Unless the L2 box contains a CKM 
phase, the imaginary part coming from RPV is zero, and 
there is no bound coming from ek] the only constraint is 
that coming from AmK- These entries have been marked 
with a dagger. 

For the B'^ — system, one noteworthy thing is that 
the bounds on the real and the imaginary parts of any 
product coupling are almost the same. (That is why we 
show only the real part: the same bound applies to the 
imaginary part too. We found it to differ by at most 
10% in a few cases.) This is, of course, no numerical 
accident. To understand this, let us analyze the origin 
of these bounds. There are two main constraints for the 
B system: Ams and sin(2/3). There will be a region, 
centered around the origin (since Am^ can be explained 
by the SM alone) of i?e(A'A') - /to(A'A') plane, where 
I A' A' I is small and the phase can be arbitrary. This is the 
SM-dominated region, where RPV creeps in to whatever 
place is left available. A conventional analysis, taking 
both SM and RPV but assuming incoherent sum of am- 
plitudes, should generate this region only. 

However, there is always scope for fully constructive 
or destructive interference. Consider a situation where 
the RPV contribution is large, so large that even af- 



TABLE II: Bounds on A'A' combinations from B° - W mix- 
ing. Rest is same as table 1, except that the entry marked 
with a dagger tree-level B° — mixing. 



ter a destructive interference with the SM amplitude, 
enough is left to saturate Ams- This RPV-dominated 
region (this is true for all NP models in general) gives us 
the bounds, and in the limit where the SM can be ne- 
glected, the bounds on i?e(A'A') are almost the same as 
on /m(A'A'). The bounds are, however, slightly different 
for the case where the RPV couplings are all real to the 
case where they can be complex; the reason is the SM 
phase in mixing, which chooses a particular direction in 
the Re{\'\') — Ini(X'X') plane not coinciding with the 
real axis. This complication is absent for the K system; 
the SM amplitude is almost real. The pattern can be un- 
derstood from Figures 2(a)-2(b). It is obvious that the 
bound is controlled by the Atob data. In Fig. 2(a) for 
the black boxes, the central cross is the region where L2 is 
important (and the interference with SM); however, the 
L4-dominated region gives the bound. Fig. 2(b) shows a 
similar scan plot for the K° — K'^ box driven by A-3iA^32, 
but without ek (this also gives the reader an idea how 
strong the ek constraint can be). However, this analysis 
holds even for the cases with only L2 boxes. 

The relative magnitude of the bounds is also easy 
to imderstand. Consider, for example, the bounds on 
A^3]^A^23 vis-a-vis A^ji^iss- The relevant box diagrams 
have the same particle content; but the first one is pro- 
portional to VibVcd A), and the second one to VtdVcb 
(~ 0{\^)). The relative suppression in A enhances the 
limit on the RPV coupling. 

It is interesting to note that in most of the cases, the 
bounds are far better than the DPBs. The exceptio n is 
the product bounds responsible for tree-level — 
and B^ — B^ mixing [9]. Apart from them, A-^^A-is has 
a comparable bound from B^ — > tt+tt" [30]. 

Though mostly of the same orders of magnitude, these 
bounds are theoretically an improvement over those ob- 
tained earlier [13-15]. We have taken into account all 
possible amplitudes (and the interference patterns play a 
nontrivial role), including the SM one, but have system- 
atically neglected the pure supersymmetric boxes coming 
from gaugino exchange. The reason is that those boxes 
decouple in the heavy squark limit, and one can always 



9 



,002 -0.0015 -0.001 -0.0005 0.0005 0.001 0.0015 0.002 
Re[(i31)(i33)] 

(a) 



0.004 p 
0.003 

0.002 - 
^ 0.001 - 

CO 

S- - 

^-0.001 - 
-0.002 
-0.003 - 
-0.00, 



3 




*Ao04 -0.003 -0.002 -0.001 0.001 0.002 0.003 0.004 
Re[(i31)(i32)] 

(b) 



FIG. 2: (a) Allowed parameter space for X'isiX'i^s. The crosses indicate the allowed parameter space switching the sin(2/3) 
constraint off, and the black boxes do the same by keeping it, along with Ams. (b) The parameter space for \i31\i32, which 
drives the — K° mixing. Only ArriK constraint is used. 



take an RPV model embedded in a minimal supersym- 
metric theory where such FCNC processes are somehow 
forbidden. It was shown in [15] that the bounds are fairly 
robust even if one takes into account such SUSY contri- 
butions. Furthermore, the QCD corrections are imple- 
mented upto NLO. We have also corrected a sign mistake 
(which, however, did not affect the result of the earlier 
papers since there was no interference to be considered) 
in the respective cfFcctive Hamiltonians. Lastly, wc; have 
also put bounds on the imaginary parts of the couplings. 
As we have stressed earlier, the role of the sin(2/3) data 
is to banish a sizable portion of the Re(X'X') — /m(A'A') 
plane, rather than constraining the absolute bounds, see 
Fig. 2(a). For the K system, ek chooses the region about 
the origin. 

If we take a baryon-number violating RPV model, 
the only combinations that can be constrained from the 
box amplitudes arc X^j^2^i2a (from B'^ — mixing) and 
^ii3'^i23 (from K'^ — if mixing). The bounds are weaker 
by about two orders of magnitude from those on their A' 
counterpart. 



V. CORRELATED CHANNELS 

Though this paper is mainly on constraints on the A' A' 
product couplings coming from K° — and — B^ 
boxes, let us also mention that such couplings are also re- 
sponsible for squark-mediated semileptonic (6 d£~^i~, 
s de+£-) and purely leptonic {B° -> £+£-, K° 
£'^£~) B and K decays (and slepton- mediated nonlep- 
tonic decays) wherever kincmatically possible. 

For the B mesons, no such leptonic mode has yet been 
observed. The corresponding upper limits on the BRs are 
of the order of 10~^ for £ = e, /i and 10~^ for the r modes. 
With 300 GeV squarks, from the bounds that one obtains 



here, a BR at most of the order of 10~ can be expected. 
Thus, we do not envisage to see such leptonic channels be- 
fore the next-generation hadronic or super e+e" B facto- 
ries. The semileptonic modes B K'-*H^£^ have been 
observed. However, the BRs are at the SM ballpark, and 
the RPV contributions are expected to be smaller by at 
least one order of magnitude. 

We have not discussed the Bg system. That will evi- 
dently warrant a detailed analysis once the hadronic ma- 
chines start running. Let us mention that the chan- 
nel Bg — > A'~'~At~ can be mediated by a tree-level 
squark exchange diagram, and with the present bound 
of BR{Bg ^ ii+jjL-) < 0.95 X 10-*^ [1], the relevant up- 
per limit on the RPV product coupling X'2i2^'2a is about 
6 X 10~^ (for 300 GeV squarks). The limit obtained from 
the Bs — Bg mixing should be at the same level after a 
few years of data. Thus, one can have an interesting sit- 
uation where there is a possibility of NP in both mixing 
and decay. 

The situation in the K meson system is better as far 
as the leptonic and semileptonic modes are concerned. 
First, note that the four- Fermi effective Hamiltonian for 
leptonic and semileptonic decays is of the form [29] 

HrPV = ^Bjklm [tjl^Pbil] [dml^PRdk] 
+ ^Bjklm \pjl'^Pm\ [dmliJ.PRdk] 
-^Cjklm M-Y^PlI^i] \dk7iJ.PLdm] + H.C., 

(37) 

where 

3 \' *\' 3 \/ *y 

Djklm - 2^ ' ^Oklm - • ^"^^1 

i=l "/rfi,. *=1 d[i^ 
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Wc take any one to be nonzero at a time. Assuming 
the leptonic phase-space to be the same, and the if + 
■K~^vV to be dominated by RPV SUSY, one obtains [10] 



Br{K+ n+vV) 
Br{K+ TfOe+z/) ~ AGprrAv:^, 



(39) 



Taking Br(Ke3) = 
n+vV) = 3.4 X 10 



= 0.0481 (lower hmit) and Br{K+ 
{la upper hmit), one obtains 



|A^^iA^2MA^yA^2^l< 3.9x10- 



(40) 



It is clear that the modes 7r2e or 7r2/i will yield a less 
severe bound, since the upper limit is weaker by orders 
of magnitude. 

The decay bounds the imaginary part 

of the product \'2ii^'2i2- The reason is that in the limit 
of CP conservation, Kl = {K^ — ii'°)/\/2, and from eq. 
(37), it is easy to sec that the real parts cancel out if 
we take the Hermitian conjugate term into account. The 
BR is 



Br{KL 



Im{X'2iiX'2i2 




(41) 



where tkl is the lifetime of K^. Putting the values, the 
bound is 



/m(A2iiA2i2) 



< 3.5 X 10" 



(42) 



This is, of course, an approximate bound, since the SM 
has not been included, and vacuum insertion is question- 
able for a t-channel squark exchange diagram. 

In the K system, the nonleptonic decays are riddled 
with theoretical uncertainties. That is why one can al- 
ways evade the bounds obtained from, say, e'/e [11]. On 
the other hand, we do not have enough nonleptonic data 
from the B factories to put comparable bounds on the 
couplings we have considered here. The only exception is 
AaiAii3, which can be bounded from the BR and the CP- 
asymmetry data of B ^ TT+TT". Other affected modes are 
listed in Table 3. However, we note that unless a product 
coupling is at least of the order of 10~^, the RPV con- 
tribution is unlikely to affect the SM amplitude. Still, a 
systematic study of NP is worthwhile; for example, the 
CP-asymmetry in _B ^ Xstt", which is supposed to yield 
sin(2/3), gives a smaller central value [31], notwithstand- 
ing the fact that the experimental errors are large: 



sin(2/?)(S ^ ifs7r°) 



0.48t°-j^ ± 0.06. 



(43) 



Note that this channel is mediated by a coupling which 
has a bound of the order of 10~^ only. 
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TABLE III: Some of the possible nonleptonic transitions me- 
diated by the RPV couplings discussed in the paper. Wc show 
only the more interesting quark-level decays in column 2, and 
some typical decays in a generic manner in column 3. The 
last column shows the UT angle that may be affected due to 
the presence of RPV. 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we have computed the bounds on the 
product of two RPV couplings of the type A'A', com- 
ing from K'^ — and — mixing, as well as from 
the sin(2/3) constraint. Though such a calculation is not 
new, we have implemented several features in the anal- 
ysis which have not been taken into account in earlier 
studies. We have considered the exact expression for the 
box amplitudes, and have taken all possible amplitudes, 
including that of the SM, into consideration. The QCD 
corrections to the amplitudes have been taken upto the 
NLO level. We have considered the possibility that the 
RPV product couplings may be complex. The analy- 
sis is done in the benchmark point mH+ = 200 GeV, 
tan/3 = 5, all sleptons degenerate at 100 GcV and all 
squarks degenerate at 300 GeV, and neglecting the pure 
MSSM contribution to the box amplitudes (by possibly 
applying to some underlying FCNC suppression princi- 
ple, like alignment of the squark mass matrices). 

It is to be observed that in some cases, our bounds are 
actually weaker than those obtained earlier by saturating 
the mass difference with RPV alone. The reason is that 
destructive interference with the SM amplitude plays a 
very crucial role in determining the bounds, particularly 
when the phase of the RPV coupling is arbitrary. There 
is an intricate interplay among different amplitudes as 
can be seen in Fig. 2(a). All in all, we consider these 
bounds to be the most conservative (or in other words 
most robust) ones. 

Some of these bounds can, however, be bettered if 
we consider the leptonic and semileptonic K decays. Of 
course, one can enhance the squark mass to a limit where 
these bounds become weaker than those obtained from 
the box (the latter is not much affected by decoupling 
the squarks) , but such extremely massive squarks arc not 
interesting probably even for the Large Hadron Collider 
(LHC). There is no such competitive bound for the B sys- 
tem, and one must wait for the hadronic B machines, or 
the super e+e~ B factories. However, some of these cou- 
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plings may affect the nonleptonic decay modes (which, 
being slcpton mediated, cannot be suppressed by decou- 
pling the squarks) of B, most important of them being 
B TT^Tr~, which has been dealt with earlier [30]. 

The Bg system is a different proposition. With the first 
data on AMg, one should try to compute all the possi- 
ble ways NP can affect this. Furthermore, one should 
also try to see whether such contributions to mixing also 
naturally lead to sizable (and comparable with the SM) 
contributions in decays like B^ J/tp(l), ^(/>, A''''M~ 
et cetera. In this respect RPV emerges as an excellent 
prospective candidate. 
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APPENDIX A: THE KINEMATIC FUNCTIONS 
COMING FROM THE BOX INTEGRALS 

The kinematic functions for the L2 boxes have been 

computed in the 't Hooft Feynman gauge, and we have 
explicitly checked that it is gauge independent. 

First, let us deal the case where there is only one type 
of charge +2/3 quark in the loop. We denote it by q. The 
bosons inside the box are a slepton (generically denoted 
as I) and one among the W boson, the charged Higgs 
boson (of mass nih) and the charged Goldstone boson 
(of mass mw in the 't Hooft- Feynman gauge). 

To start with, let us define some shorthand notations: 
Q = ml, L = m?, W = m^, if = ml (Al) 
The kinematic factors are 



" {L-W){L-Q)^W-Q)^ 

[Q{L - W){2LW - Q{L + W)) logQ + L'^{W - Qf logi 
+ (L - Q){Q{W - Q){L -W)- W^L - Q) logW}], 
a _ [W{L - Qf log W+{L- W){Q^ - LW) log Q - {W - Q){{L - W){L - Q) + L{W - Q) log L}] 
~ {L-W){W -Qf{L-Qf 

■A^2 = A%{W ^ B). (A2) 

The Hamiltonian is slightly modified if we have two different up-type quarks q\ and 92. Using the shorthand Qxifi'i) = 
2 ' one has 

'Sfa = ;^ [Qi J^i log gi + Q2G1 log g2 + LJ^ log L - WK^ log W] , 

= --^[FilogQi + GilogQ2 + JilogL + KilogVF], 
Bf2 = B%{W ^ H). (A3) 



Here 



D = {Qi-Q2){L-Qi){L-Q2){W -Q{){W -Q2){L-W), 
Fi = Q^{L-Q2){W -Q2){L-W), 

Gi = -Fi(Qi ^02), 
Ji = {Qi-Q2){W-Q^){W-Q2)L, 

Ki = -Ji{L ^ W). (A4) 



We have checked that in the limit Qi = Q2, the B APPENDIX B: THE FOUR- A' BOXES 

functions reduce to corresponding ^s. 



Apart from the boxes that we have discussed, there 
may be boxes whose amplitudes are proportional to the 



12 



4A 


Internal 


Operator 


combination 


propagators 






same u, different £ 


Ui 




same d, different v 


Oi 




same u, different £ 






same d, different f 


di 


{ilj)(klj){iqj) (kqj) 


same d, different u 


Ul 




same d, different u 


r\ 
U\ 


{ijq)(ikq){ijl) (ifcl) 


same different u 


Ul 




same z/, different d 


Oi 




same £, different u 


Oi 




same i>, different d 


Oi 


{ilj)(ilk)(iqj) [iqk) 


same z/, different d 


Oi 




same i>, different d 


Oi 


{ijq){lkq)[ikl) (Ijl) 


different u and £ 


Ul 




different d and 


Ul 




different f and u 


di 




different u and d 


6i 


{ilk){llj){iqjr{lqkr 


different d and /> 


Oi 




different v and d 


Oi 



TABLE IV: Different 4A' combinations, and the correspond- 
ing propagators. 



product of four different A' couplings. They are not 

important if wc assume a hierarchical couphng scheme 
where one and only one product is numerically signifi- 
cant, but may appear in a scheme with a given texture 
at a high scale which results in a number of nonzero cou- 
plings at the weak scale. Some of these couplings have 



been shown in, e.g., [14]. 

The effective Hamiltonian is always proportional to, 
apart from the four A' couplings, a factor of l/647r^ times 
a momentum-space integral of the form 



F(/) 



+ m?) 



(Bl) 



where the denominator indicates the four propagators 
inside the box. To obtain a quick numerical estimate, 
we may put all SM fermions except the top quark to be 
massless, all sleptons (neutral and charged) to be degen- 
erate (say at 100 GcV) and all squarks to be degenerate 
too (say at 300 GeV). The operator is either Oi or Oi, 
see eq. (14). If the sneutrinos are not their own antipar- 
ticles, there cannot be any operators of the form O4 or 
O5. 

In the following table, we show all possible combi- 
nations of the 4A' product, with the propagators. Ex- 
panded, they look like, e.g., for same up-type quarks but 
different sleptons: 



A' A' A'* A'* 

ijq kjq ij^ kjl ^ 

647r2 



OiF(p2) 



(B2) 



where F(p^) is to be evaluated with two quarks and two 
different sleptons. 

Note that if these combinations arc present, bounds 
coming from L2 boxes alone have to be reevaluated. 
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